Self-Calibrated Cluster Counts as a Probe of Primordial Non-Gaussianity 
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We show that the ability to probe primordial non-Gaussianity with cluster counts is drastically 
improved by adding the excess variance of counts which contains information on the clustering. The 
conflicting dependences of changing the mass threshold and including primordial non-Gaussianity 
on the mass function and biasing indicate that the self-calibrated cluster counts well break the 
degeneracy between primordial non-Gaussianity and the observable-mass relation. Based on the 
Fisher matrix analysis, we show that the count variance improves constraints on /nl by more 
than an order of magnitude. It exhibits little degeneracy with dark energy equation of state. We 
forecast that upcoming Hyper Suprime-cam cluster surveys and Dark Energy Survey will constrain 
primordial non-Gaussianity at the level ct(/nl) ~ 8, which is competitive with forecasted constraints 
from next-generation cosmic microwave background experiments. 

PACS numbers: 95.36.+X, 98.65.Cw, 98.80.-k 



The measurement of departures from Gaussianity of 
the initial perturbations provides a unique opportunity 
to probe the early universe [!}. While the standard sin- 
gle field slow-roll inflation models predict primordial per- 
turbations very close to Gaussian, some models such as 
multi-field models and the curvaton model can produce 
the level of non-Gaussianity high enough to be detected 
in ongoing or future surveys. Thus specific forms of pri- 
mordial non-Gaussianity contain valuable information on 
how the initial density fluctuations are generated. 

Observationally, primordial non-Gaussianity has 
mainly been studied using the temperature fluctuation 
of the Cosmic Microwave Background (CMB). Recently 
it has attracted considerable attention given a possible 
detection of non-Gaussianity by Yadav and Wandelt [2| • 
However, the detection of non-Gaussianity in the CMB 
is somewhat controversial in the sense that independent 
analyses yield slightly different results Q , suggesting the 
importance of other observational probes independent 
of the CMB. Another powerful probe of primordial 
non-Gaussianity is provided by the large-scale structure 
of the universe. In particular, non-Gaussianity induces 
a scale-dependent halo bias 0, H, 0, > and thus by 
studying large-scale distributions of astronomical objects 
one can obtain tight constraints that are competitive 
with the CMB. Constraints from the large-scale struc- 
ture are also important given that non-Gaussianity can 
be scale-dependent such that deviations from Gaussian 
are larger at smaller scales @. 

Primordial non-Gaussianity is also sensitive to the 
abundance of massive clusters and its redshift evolution 
0, flol | ■ An advantage of using massive clusters is its 
reasonable one-to-one correspondence with dark halos, 
which suggests that halo assembly bias (e.g., Q) is less 
important. A challenge here is how to calibrate clus- 
ter masses; since the cluster mass is not directly observ- 
able, one has to resort to well-calibrated correlations be- 
tween cluster masses and observable quantities such as 



luminosities, temperatures, and the numbers of member 
galaxies in order to infer cluster masses. The observable- 
mass relations always involve uncertainties, suggesting 
that the change of cluster abundances by primordial non- 
Gaussianity may be compensated by modifying the re- 
lation between observables and masses. Therefore con- 
straints from cluster counts rely on how well we can cal- 
ibrate such observable-mass relations. 

In this Letter, we point out that clustering informa- 
tion breaks the degeneracy and allows us to determine 
primordial non-Gaussianity surprisingly well with clus- 
ter counts. This is because the clustering bias for mas- 
sive clusters is quite sensitive to both cluster masses and 
primordial non-Gaussianity, and more importantly, be- 
cause the cluster abundance and biasing show conflicting 
dependences on these. Such self-calibrated cluster count 
technique has been discussed extensively in the context 
of accurate dark energy probes [ill, 112, |l3, LLJ] , but its use 
for primordial non-Gaussianity has not been explored. 

Here we quantify non-Gaussianity of the local 
form using the standard parametrization, $ = (f> + 
/nl (<t> 2 — (4> 2 ))i where <J> is the curvature perturbation 
and <f> is an auxiliary random-Gaussian field. The pa- 
rameter /nl > (< 0) indicates that the initial den- 
sity field is positively (negatively) skewed. The cur- 
rent level of constraints on primordial non-Gaussianity 
is 1/ 



NL 



< 



0(100) [2, S We adopt a non-Gaussian 
correction factor of the cluster mass function based on 
the Edgeworth expansion [§]: 
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where v = S c /a, 5 C w 1.68 is the critical linear overden- 
sity, o~ — cr(M, z) is the linear fluctuation on the mass 
scale of M which we compute using the transfer function 
T{k) presented by Eisenstein and Hu [15| ignoring the 
baryon wiggle. We adopt models of Warren et al. [l6| 
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for the mass function in the Gaussian case, dnc/dM. 
The skewness S3 is related to /nl as 

aS3 = ^ > ^ia(fc 1 )W(M,fc 1 )A2(fc 1 ) 
17 Jo ft i 
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where A| = k 3 P 4> {k)/27T 2 is the power spectrum of the 
curvature perturbation, a = [2Z)(z)T(fc)/3f2M]( c ^/^o) 2 , 
D(z) is the linear growth rate normalized to (1 + z)^ 1 in 
the matter-dominant era, and k 2 = k 2 + k 2 + 2/j,kik2- For 
the window function W(M, k) we adopt the real space 
top-hat filter. In practice we use the following fitting 
formula for a S3: 
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^b(l + v2/i/57m)] is so-called the shape pa- 



17] . This fitting formula should be accurate at 
a few percent level in the mass scale range 10 7 /i _1 M Q < 
M < IO^/i-^q. 

The non-Gaussian correction of the halo bias is com- 
puted as m 



Ab{M,z,k) = 



a 



]' d ( dn/dM \ 
5 C dv 



^ G ^ s J -- [dn G /dM) 



(4) 

The halo bias in the Gaussian case, be, is assumed to be 
the form presented by Sheth and Tormen [l8j ]. 

Fig.Q]illustrates the reason why the clustering informa- 
tion is so important. As shown in the Figure, including 
positive /nl increases the number of clusters above some 
mass threshold M t h- This increment can be compen- 
sated by raising M t h- However, these two models with 
the same numbers of clusters result in quite different halo 
biases because both raising M t h and /nl increase the bi- 
asing. Thus by including clustering information we can 
strongly break the degeneracy between M t h and /nl, and 
can obtain tight constraints on /nl- 

We now forecast constraints on /nl from future cluster 
surveys. We include the clustering information using a 
count-in-cell analysis. Specifically we approximate the 
Fisher matrix as 12. 1 1 311 



0-2 (a)' 



(5) 

where <r p represents the prior information on each pa- 
rameter, and the covariance matrix is given by C = 




M [h^MJ 



FIG. 1: The illustration of self-calibrating cluster counts to 
probe primordial non-Gaussianity /nl- Here the redshift is 
z = 1. In cluster surveys, we basically obtain the number 
of clusters above some mass threshold Mth- The plot in- 
dicates that the increase of the mass function dn/dM due 
to positive /nl can be compensated by increasing the Mth- 
However, while the two models predict the same number of 
clusters, the corresponding halo bias b(M) (here we adopted 
k = 0.02/iMpc" 1 ) are quite different because both raising 
Mth and /nl increase 6(M) (see arrows). Although in this 
plot the scatter in the observable-mass relation is ignored for 
simplicity, it will be included in the Fisher matrix analysis. 



S + diag(m). The number count m and its variance S 
are computed as 
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where the subscript i run over redshift, mass, and an- 
gular bins. The power spectrum is described by P(k), 
and the fc-space window function by Wi(k). Since the 
off-diagonal elements of S are small in our case, here 
we consider only the diagonal elements. The function 
p(M b s \M) models the accuracy of the cluster mass de- 
termination from observables. Following [l3j |. we assume 
the log-normal distribution for p{M \ is \M) 1 with the me- 
dian of In M + In Mbias and the scatter of a\ n m , and re- 
gard cin m and In Mbias (which corresponds to Mth in 
Fig. [T]) as nuisance parameters. Note that the first term 
of the Fisher matrix (Eqn. [5]) represents the informa- 
tion from number counts, whereas the second term the 
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information from the variance of the counts which con- 
tain clustering (biasing) information. Using the Fisher 
matrix, one can estimate a marginalized error on each 
parameter as a (a) = y/ (F -1 )^. 

We calculate the Fisher matrix in 10-dimensional pa- 
rameter space; 6 standard cosmological parameters in- 
cluding dark energy equation of state (the matter den- 
sity SIm'i 2 , the baryon density flbh 2 , the power spec- 
trum tilt n s , the normalization of the power spectrum 
<^c (3> the dark energy density JIde, and dark en- 
ergy equation of state w), 1 parameter representing 
primordial non-Gaussianity (/nl), an d 3 parameters 
from the observable-mass relation, o~\ n M and lnMbi as = 
InMbias.o + 7ln(l + z). The Five- Year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP5) result for ACDM 
[H, (fl M h 2 , n b h 2 , n s , S c , n UE , tu)=(0.133, 0.0227, 
0.963, 4.61 x 10~ 5 , 0.742, -1), is adopted as our fidu- 
cial cosmological model. We add conservative priors to 
the first 4 parameters, a v (VLMh 2 ) = 0.006, a p (flbh 2 ) — 
0.0006, a p (n s ) = 0.015, and a p (6 c ) = 10~ 6 ; these arc the 
level of accuracies which has already been achieved by 
WMAP5. In addition, our fiducial model has /nl = 0, 
cinAf = 0.25, In M b ias,o = 0, and 7 = 0. 

For illustrative purposes, we consider the following 
three upcoming surveys; Hyper Suprime-Cam on Sub- 
aru telescope (HSC; since the design of the HSC cluster 
survey is still tentative, we consider both 1000 deg 2 and 
2000 deg 2 ), Dark Energy Survey (DES; 5000 deg 2 ) [Hj], 
and Large Synoptic Survey Telescope (LSST; 20000 deg 2 ) 
[22] ] . We adopt a simplified assumption that these optical 
surveys will find clusters with M Q bs > 10 13 ' 7 h~ 1 M (3 out 
to z max = 1.4, 1.0, and 1.7, respectively. For the count- 
in-cell analysis, we use the cell size of 20 deg 2 and the 
redshift interval Az = 0.1. Three mass bins with spacing 
of AlogAf bs = 0.5 are also adopted. 

In Fig. [2j we show marginalized constraints on /nl 
and the correlations with the parameters a\ n m , expected 
for the 2000 deg 2 HSC cluster survey. As expected, con- 
straints are drastically improved by combining the num- 
ber counts with the variance which includes the clus- 
tering information. This is partly because constraints 
from number counts and variance show different degen- 
eracy directions, suggesting that both the number counts 
and clustering are essential for accurate determinations 
of /nl- 

Table U summarizes forecasted constraints on various 
cosmological parameters. For all the upcoming surveys, 
the count variance drastically enhances the ability to 
probe primordial non-Gaussianity, by more than an or- 
der of magnitude improvement in /nl compared with 
the number counts alone. Predicted marginalized errors 
of <t(/nl) ~ 8 for HSC and DES and - 2 for LSST are 
competitive with constraints from next-generation CMB 
experiments (e.g. . |23| ) and galaxy power spectrum mea- 
surements (e.g., [4j, |24|). The variance helps to regulate 
the observable-mass relation, improving an accuracy of 
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FIG. 2: Expected marginalized constraints in the /nl-ciuM 
plane from the upcoming 2000 deg 2 HSC cluster survey. 
WMAP5 cosmology is assumed as a fiducial model. Con- 
tours indicate the 68% confidence regions from the number 
counts (dotted), the variance of the counts which contains the 
clustering information (dash-dotted) , and the combination of 
the number counts and variance (solid). 



InMbias.o and a\ n M measurements by a factor of two or 
more. Measurements of dark energy equation of state are 
improved as well, which is consistent with earlier work. 
We find that w and /nl are not correlated very much, in- 
dicating that we can well determine these two parameters 
simultaneously using self-calibrated cluster counts. Fig. [3] 
shows contours of a (/nl) as a function of the survey area 
and the maximum redshift. The expected constraints on 
/nl is a steep function of the maximum redshift even at 
z > 1, which indicate the importance of the deep surveys 
to detect clusters out to z > 1. 

We have shown that adding clustering information 
from the count variance drastically improves measure- 
ments of primordial non-Gaussianity with cluster counts. 
Although the calibration of cluster masses limits the 
use of cluster counts as a cosmological probe, the self- 
calibration technique allows us to determine both the 
observable-mass relation and /nl simultaneously. The 
significant effect of the count variance comes from the 
conflicting dependences of the mass threshold and /nl on 
the cluster mass function and biasing parameter (Fig.[T]). 
Allowing dark energy equation of state to vary does not 
degrade /nl measurements very much. Resulting fore- 
casted constraints on /nl, c(/nl) ~ 8 for HSC and DES 
and ~ 2 for LSST, suggest that cluster counts can be- 
come a competitive probe compared to the CMB or the 
large-scale galaxy power spectrum. 

We have here made a number of simplified assump- 
tions. For instance, it is important to check how the pos- 
sible redshift evolution of the observable-mass relation 
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Survey 


Counts 


Variance 


Counts + Variance 


o-(fiDE) cr(w) <t(/nl) 


ct(Qt)e) ar{w) <t(/nl) 


a(fl OE ) (t(w) <t(/nl) 


HSC1 
HSC2 
DES 
LSST 


0.030 0.151 240.4 
0.023 0.108 188.2 
0.032 0.081 210.6 
0.009 0.037 106.1 


0.012 0.189 37.8 
0.011 0.142 28.0 
0.011 0.102 35.3 
0.006 0.051 6.9 


0.010 0.103 8.1 
0.009 0.074 6.2 
0.009 0.055 7.9 
0.005 0.024 1.9 



TABLE I: Marginalized constraints on cosmological parameters estimated from the Fisher matrix analysis using the number 
counts and/or the variance of counts (clustering). WMAP5 cosmology is assumed as a fiducial model. Constraints in four 
future survey parameters, HSC1 (1000 deg 2 , z max = 1.4), HSC2 (2000 deg 2 , z max = 1.4), DES (5000 deg 2 , z max = 1.0), and 
LSST (20000 deg 2 , z max = 1.7), are presented. 
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FIG. 3: Expected marginalized constraints on /nl as a func- 
tion of two survey parameters, the survey area and the max- 
imum redshift z max . Contours are drawn per 0.25 dex. Four 
survey parameters considered in this paper are indicated by 
filled squares. 



affects our results [12|, |13[ . The impact of other system- 
atics, such as cluster photometric redshifts [14| and the 
effect of halo assembly bias [25j], should be addressed. 
On the other hand, we used only the count variance of 
each cell as the clustering information. Since the effect of 
primordial non-Gaussianity is more significant at larger 
scales, including the count covariance (or including the 
full clustering information with the power spectrum (ll| ) 
may improve the constraints further. We leave such more 
comprehensive treatments for future work. 

This work was supported by Department of Energy 
contract DE-AC02-76SF00515. 
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